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A method to construct exact general relativistic thick disks that is a simple generalization of the 
"displace, cut and reflect" method commonly used in Newtonian, as well as, in Einstein theory of 
gravitation is presented. This generalization consists in the addition of a new step in the above 
mentioned method. The new method can be pictured as a "displace, cut, fill and reflect" method. 
In the Newtonian case, the method is illustrated in some detail with the Kuzmin-Toomre disk. 
We obtain a thick disk with acceptable physical properties. In the relativistic case two solutions 
of the Weyl equations, the Weyl gamma metric (also known as Zipoy-Voorhees metric) and the 
Chazy-Curzon metric are used to construct thick disks. Also the Schwarzschild metric in isotropic 
coordinates is employed to construct another family of thick disks. In all the considered cases we 
have non trivial ranges of the involved parameter that yield thick disks in which all the energy 
conditions are satisfied. 
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I. INTRODUCTION 

Exact solutions of the Einstein equations are associ- 
ated to highly idealized physical systems that have some 
exceptional geometrical properties. In some cases with a 
simple exact solution one can capture a significant part 
of the physical properties of non trivial systems. Also, in 
nonlinear theories like general relativity and fluid dynam- 
ics the exact solutions play an important role in numeri- 
cal analysis. These solutions can be used to test numer- 
ical codes and its outcomes. Also they can be employed 
as initial conditions to describe more realistic situations, 
e.g., a static solution can be used as part of the initial 
conditions for a full dynamical simulation. 

Since the natural shape of an isolated self-gravitating 
fluid is axially symmetric, the solutions of Einstein's field 
equations with this symmetry play a particularly impor- 
tant role in the astrophysical applications of general rel- 
ativity. In particular disk like configurations of matter 
are of great interest, since they can be used as models 
of galaxies or accretion disks. Also these disks can be 
used as starting point to represent more realistic models 
in which the bulge and halo of the galaxy are considered. 

Solutions for static thin disks without radial pressure 
were first studied by Bonnor and Sackfield P| , and Mor- 
gan and Morgan , and with radial pressure by Morgan 
and Morgan Q. The first solution represent disks made 
of pressureless dust whereas the second disks with az- 
imuthal pressure but without radial pressure. And the 
third a disk made of an anisotropic fluid with nonzero ra- 
dial pressure. The Bonnor-Sackfield disk has a singular 
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rim. These disks are finite. 

Several classes of exact solutions of the Einstein field 
equations corresponding to static 4c i& 7, 8, 9, 10, l3| 
and stationary 0, 0, Q thin disks have been obtained 
by different authors, with or without radial pressure. 
Thin disks with radial tension ITsll . magnetic fields 
and magnetic and electric fields 17] have been also stud- 
ied. The non linear superposition of a disk and a black 
hole was first considered by Lemos and Letelier ■ This 
solution and its generalizations has been studied in some 
detail in ^ia,^2L22»^,lil5j. Recently the sta- 
bility of circular orbits of particles moving around black 
holes surrounded by axially symmetric structures have 
been considered in |26j |. For a recent survey on relativis- 
tic gravitating thin disks, see j27j. 

Except for the pressureless disks all the other disks 
have as source matter with azimuthal pressure (tension) 
different from the radial pressure (tension). However, in 
some cases these disks can be interpreted as the superpo- 
sition of two counter-rotating perfect fluids. A detailed 
study of the counter-rotating model for the case of static 
thin disks is presented in |28|| . Recently, more realistic 
models of thin disks and thin disks with halos made of 
perfect fluids were considered in [29^ . 

In all the disks mentioned above an inverse style 
method was used to solve the Einstein equations. The 
metric representing the disk is guessed and then it is used 
to compute the source (energy- momentum tensor). This 
method was named by Synge as the g-method j30l | in 
contrast to the t-niethod or direct method in which the 
source is given and the Einstein equations are solved. The 
t-method, has been used to g enerate disks by the Jena 
group m m m Hi m IM 113. Essentially, they are 
obtained by solving a Riemann-Hilbert problem. These 
solutions are highly non trivial, but they deserve special 
attention because of their clear physical meaning. 
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In the solutions obtained by the g-method the well 
known "displace, cut and reflect" method is used. The 
idea of the method is simple. Given a solution of the 
vacuum Einstein equations, a cut is make above all sin- 
gularities or sources. The identification of this solution 
with its mirror image yields relativistic models of disks. 
In general, these disks are of infinite extension and finite 
mass. 

The aim of this paper is to consider disks beyond the 
thin disk limit to add a new degree of reality to these 
geometric models of galaxies. Even though in first ap- 
proximation the galactic disks can be considered to be 
very thin , e.g., in our Galaxy the radius of the disk is 
lOkpc and its thickness is Ikpc. In a more realistic model 
the thickness of the disk need to be considered. Also it is 
well known in fluid mechanics that the addition of a new 
dimension can make dramatic changes in the dynamics 
of the fluid. In principle, this new dimension will also 
change the dynamical properties of the the disk source, 
e.g., its stability. 

In this paper we generalize the "displace, cut and re- 
flect" method in order to obtain thick disks models from 
vacuum solutions of Einstein equations. We shall replace 
the surface of discontinuity of the metric derivatives with 
a "thick" shell in such a way that the matter content of 
the disk will be described by continuous functions with 
continuous first derivatives. This generalization can be 
named "displace, cut, fill and reflect" method. The disks 
obtained with this method, in general, will be of infinite 
extension and finite mass. Also as in the case of thin 
disks, the matter that form the thick disks will not obey 
simple equations of state and in some regions of the disk 
the pressure can change sign given rise to tensions. Al- 
though the energy condition will be fulfilled. The models 
of thick disks presented can be considered aS generaliza- 
tions of models of thin disks studied in references 
and m. 

The article is divided as follows. In Sec. II we present, 
in some detail, the main idea of the "displace, cut, fill 
and refiect" method in Newtonian gravity. The method 
is then applied, in Sec. Ill, to construct relativistic thick 
disks in Weyl coordinates. We also study the general ex- 
pression for the energy-momentum tensor of the disks. 
The method is illustrated by taking two simple Weyl so- 
lutions that lead to thick disks in agreement with all the 
energy-conditions. In Sec. IV we apply the method to 
the Schwarzschild solution in isotropic cylindrical coordi- 
nates. The disk obtained also satisfy all the energy con- 
ditions, this disk have equal azimuthal and radial pres- 
sures and different vertical pressure. Finally, in Sec. V, 
we summarize our main results. 



II. NEWTONIAN THICK DISKS 

The Newtonian gravitational potential of a thin disk 
can be obtained by a simple procedure, the "displace, cut 
and reflect method" , due to Kuzmin [s^ . The method 



can be divided in the following steps: First, choose a 
surface that divides the usual space in two parts: one 
with no singularities or sources an the other with the 
sources. Second, disregard the part of the space with 
singularities. Third, use the surface to make an inversion 
of the nonsingular part of the space. The result will be 
a space with a singularity that is a delta function with 
support on z = 0. This procedure is depicted in Fig. ^ 

In order to obtain a thick disk we need to modify the 
above procedure. We essentially need to replace the sur- 
face of discontinuity with a "thick" shell in such a way 
that the matter content of the disk be described by con- 
tinuous functions. Now method has an additional step 
and can be named "displace, cut, fill and reflect". Af- 
ter we disregard the part of the space with singularities, 
we put a thick shell below the surface. Then we use the 
bottom surface of the shell to make the inversion. The 
procedure is illustrated in Fig. |21 

Mathematically the method is equivalent to make the 
transformation z — > h{z) + b, where 6 is a constant and 
h{z) an even function of z. In Newtonian gravity, the 
potential ^(r, z) is a solution of the Laplace Equation 

= + ^ + = 0, (1) 
r 

where (r, ip, z) are the usual cylindrical coordinates. Af- 
ter we make the transformation z — > h{z) -\-b^ the above 
equation leads to 

= h"^,h + [{h'f - 1]$,,,^, (2) 

where primes indicate differentiation with respect to z. 

For the case of thin disks we take h{z) — \z\. Note 
that d^\z\ = 26'(z)-l and d^e{z) = 5[z) where 6* (z) is the 
Heaveside function and 5{z) the usual Dirac distribution. 
By using Q , the Poisson equation leads to a mass density 
given by 

2TTGp{r,z)^'^,h5{z). (3) 

We have a surface distribution of matter located in the 
plane z = 0. In FigO with dashed lines, we plot h{z) an 
its derivatives. 

For the case of thick disks the function h{z) must be 
selected in such a way that $ and its flrst derivatives be 
continuous across the plane z — Q. Let us take a function 
h{z) deflned as 

{z — a/2 , z> 
zV2a , -a < z < a, (4) 
— z — a/2 , z < —a. 

Hence, by taking the function h{z) above deflned we can 
generate disks of thickness 2a located in the region — a < 
z < a. In Fig. Owe plot h{z) an its derivatives (solid 
lines). 

When |z| > a the fimction h{z) is a linear function of z 
such that h'{z) — 1. Hence, its second derivative is zero. 
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FIG. 1: Construction of a thin disk by the "displace, cut and reflect" method from the gravitational field of a mass point. In 
(a) the space with a singularity is displaced and cut by a plane (dashed line). In (b) the part with singularities is disregarded 
and the upper part is reflected on the plane. 




FIG. 2: Construction of a thick disk by the "displace, cut, fill and reflect" method. In (a), after disregard the part with 
singularities, we put a thick shell below the plane. In (b), we refiect the resultant configuration on the bottom surface of the 
shell. 



Then the mass density vanish outside the disk. Since the 
first derivative is continuous at |z| = a (see Fig. and 
the second derivative piecewise constant we have that the 
mass density, p, will be well defined inside the disk, 



AnGa'^pir, z) = a$,,i + {z^ - a^)$. 



hh 



(5) 



for \z\ < a and p = for \z\ > a. 

As a simple example we can consider the usual poten- 
tial for a mass point, written in cylindrical coordinates 
as 



Gm 



(6) 



By doing the transformation z —> h{z) + b in the previous 
potential, we obtain the mass densities 



p(r, z) = 



b6{z) 



27r(r2 + 62)3/2 



(7) 



for a thin disk (Kuzmin-Toomre disk |33,|39j), and 



p{r, z) 



■m[3z'^ + 2a{b-a)] 
87ra2R3 

3m(a2-z2)(^2^2Q6)2 
167ra4R5 



(8) 



where R"^ = r'^ + {h + b)^, for a thick disk. When b> a, 
the mass density will be positive everywhere. 

The function h{z) presented in Q is the simplest func- 
tion that have the desired properties. The part of the 
function in the domain |z| < a can be changed by super- 
position of even functions of z. In [zj = a this new func- 
tion need to be matched continuously with linear func- 
tions of z such that h'{z) — 1, like in 10}. In this case 
the mass density can depend on the variable z in a more 
general way. 

It is instructive to obtain the surface thin disk density 
(say a ) associated to O, 



(9) 



27r(r2 + 62)3/2 
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FIG. 3; The function h{z) and its derivatives, for thin disks (dashed line) and thick disks (solid line). In (a) we plot h{z), in 
(b) h'{z) and in (c) h"{z). 



as a limit of the thick disk volume density (jSJ. To ob- 
tain the surface density of the thin disk from the volume 
density Q) we first do 



E — 2ap{r, z). 



(10) 



Now to perform the thin disk limit we put z = and 
a — (3^ (a and f3 artitrary constants) in E and we take 
the limit lim^^o 5^. This limit is just ©• 



III. RELATIVISTIC THICK DISKS IN WEYL 
COORDINATES 



When the matter is absent, the metric for a static ax- 
ially symmetric spacetime can be cast without loosing 
generality as 



"dt' + e-^*[r^ V + e'^idr' + dz% (11) 



where $ and A are functions of r and z only. The ranges 
of the coordinates {ip, r, z) are the usual for cylindrical 
coordinates (Weyl coordinates) and —oo<t<oo. The 
Einstein vacuum equations for this metric yield the Weyl 
equations |40ll4l| . 



r 

A 7 = 2r<i> 



(12a) 
(12b) 
(12c) 



From a solution of Einstein vacuum equations corre- 
sponding to a Weyl metric (|ll|l . we can construct a thick 
disk model by means of the "displace, cut, fill and reflect 
method" using the transformation z — > h{z) + h, with 
h{z) given by Q. The energy-momentum tensor of the 



disk can be computed using the Einstein equations in the 
matter, written as 



Tab — Rah — -^gabR, 



(13) 



in units such that c = SttG = 1 . 

By using the Einstein equations H12|l . the nonzero com- 
ponents of are: 



t; = 



rpZ 





A) 




- 








A) 


a2 


- 




A) 


a2 





e2(*- 


A) 


a2 



(14a) 



[aA^h + {z' - a')($i + A,wO],(14b) 



(z^ ~ a^)<^\ , 



(14c) 
(14d) 



valid for the region < z < a. Outside this region we 
have Tl^ = 0. 

For the above expressions we can see that the radial 
stress is negative (we have radial tension). On the 
other hand, since — — TJ^^ we have vertical pressure. 
Defining the orthonormal tetrad {V'', X'', , Z''}, where 



X" 



e-* (1,0,0,0) , 



— (0,1,0,0) , 
r 



(0,0,1,0) 



(0,0,0,1) 



(15a) 

(15b) 
(15c) 
(15d) 
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we can be cast the energy-momentum tensor in its canon- 
ical form 

Tab = eKH + p^XaXb + PrYaYb + PzZaZb ■ (16) 

Here e — —T^ is the energy density, = is the 
azimuthal stress, Pr = T,r = —T^ is the radial tension 
and pz — Tl is the vertical pressure. 

From 114|) we get the "effective Newtonian" density, 

p = e -I- + + = e + 

r,p2(0-A) 

P - , + {z" - a2)<i>.,,.]. (17) 

The strong energy condition requires that p > 0, whereas 
the weak energy condition imposes the condition e > 0. 
The dominant energy condition is equivalent to the re- 
quirement 1^1 < 1, 1^1 < 1 and 1^1 < 1, see for in- 
stance [i^ . 

One can obtain the quantities associated to a general 
relativistic thin disk from the corresponding quantities 
associated to the thick disk using the same limit proce- 
dure described in the Newtonian disk case. 
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A. Thick Disks from the Chazy-Curzon Metric 



As a first example we apply the "displace, cut, fill and 
refiect" method to obtain thick disks using the Chazy- 
Curzon solution ji^, 0| , written in Weyl coordinates as 



A 



m 

R' 



2R4 



(18a) 
(18b) 



where 



(/i -I- 6)^, m and 6 are positive constants 



and h{z) is given by Q). 

Now we rescale the variables and the parameters in 
terms of the disk thickness, a. We r — af^ z ~ az, 
R = aR, b = ab and m = am. From (|14|l and (|17|l . we 
obtain 



=2(*-A) . 



2R5 




2f<E>- 

m 


-A) 


4R8 






-A) 


4R8 




^2g2(*- 


-A) 


4R6 




^2p2(*- 


-A) 



4R6 



2(3^2 + 2{b - 1))r2 + 3(1 - 52)(z2 + 2bf 

4(3z2 + 2(b - 1))(R3 - mf2)R2 + [6(r3 _ 2mr^) + mR2)(l - + 2bf 

4f2[(3i2 + 2(6 - 1))R2 + 3(1 - z2)(f2 + 2bf] + (i^ - i){z^ + 2fe)2R2 

(z2_ l)(i2 + 26)2 

"(l-z2)(z2-|-26)2 

r 



(19a) 

(19b) 
(19c) 

(19d) 

(19e) 



where p = a^p, e = a^e and pi = a^Pi- 

Equation H19a|) and the condition 5 > 1 imply p > 0. 
Then, when 6 > 1 the strong energy condition is satisfied. 
From l(T9K)l we conclude that e > whenever R"^ > 2mf ^ 
and & > 1. From the definition of R we conclude that in 
order to have e > everywhere, we need 



< m < 



^/36 



and b > 1. 

The behavior of the densities is better illustrated 
graphically. In order to have a disk in agreement with 
the weak and strong energy conditions, we take m = 1 
and b — 2. In Fig. 0] we plot the effective Newtonian 
density and the energy density in units of , p and e, as 



functions of f and z. We can see that the densities have 
a maximum at the center of the z = plane. Then the 
densities decrease monotonously as r increase and also 
the densities decrease for z — > ±a. We see that p and e 
have similar magnitudes. 

The azimuthal stress, as we can see from (I19cf) . is nega- 
tive at the center of the disk (r = 0), whereas it is positive 
for larges values of r. The boundary between the region 
of negative stress and positive stress is the surface 



4j^8g2(A-*) 



(20) 



with p^ given by (|19c() . 

The behavior of the stresses is also better illustrated 
graphically. In Fig. [S] we plot the dimensionless az- 
imuthal stress p^p and vertical pressure Pz, as functions 
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FIG. 4: For a thick disk obtained from the Chazy-Curzon solution with rh — 1 and 6 = 2 we plot (a) the dimensionless 
Newtonian density p and (b) the dimensionless energy density e, as functions of r and z. 




FIG. 5: For a thick disk obtained from the Chazy-Curzon solution with m — 1 and 6 = 2 we plot (a) the dimensionless 
azimuthal pressure and (b) the dimensionless vertical pressure pz- 



of f and z. The radial pressure is given by Pr — —Pz- 
Again we take rh — 1 and b — 2. The azimuthal stress is 
negative at the central region of the disk, then increases 
to have a positive maximum on the z — plane, for a 
value of r « 1.5a. Finally, it decreases monotonously for 
increasing r and also for z ±a. The behavior of the 
vertical pressure is like of the densities, with a maximum 
at the center of the z = plane and then monotonously 
decreasing for increasing r. Also, Pz = for z = ±a. 
From Figs. 21 and [S] we can also see that the magnitude 
of the stresses is about a tenth of the magnitude of the 
densities. We have 



Pr 




El 


e 




e 



< 0.1, 



< 0.1, 



and so the disks are also in agreement with the dominant 
energy condition. Thin disks based on the Chazy-Curzon 
metric were studied in Ref. Hi- 



B. Thick Disks from the Zipoy-Voorhees Metric 

As a second example we take the Weyl gamma metric 
also known as Zipoy-Voorhees solution (45t 1461] , that in 
Weyl coordinates can be cast as ^3 



A 



m , 
— In 

2k 



Ri 



R, - 2k 



m 
2fc2 



In 



Ri + R2 + 2k 



(Ri+R2)2-4A:2 



4R1R2 



(21a) 
(21b) 



where Rj = r'^ + {h + b + kf 



R-2 = 



+ {h + b-ky 



m, b and k are positive constants and h{z) is given by 
Q . When k = m this solution leads to the Schwarzschild 
metric and when fc to the Chazy-Curzon solution of 
the previous section. 

By using (|14() and H17() we obtain, in terms of the di- 
mensionless variables used in the previous section. 
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Pr 



1 - -TT- 



|2(Ri - R2)R?Ri + (1 - z^) [z^ + 2&)(R? - Ri) - 2fc(Rf + R^) | , 



(22a) 



2mf2(Ri +R2) 



R2R2[(Ri +R2)2 -4fc2 



m2e2(*-A)(l-z2)(j^^ _ 1^2) 



4fc2R4R4 



(Ri-R2)RfR2-2P(Rf-R^) ,(22b) 



2mpf'^{Ri+R2) m2g2(*-A)(~2_ ;l)(Ri -R2) 



R2R2[(Ri +R2)2 -4fc2] 



4fc2R4j^4 



(Ri - R2)R?R2 - 2f^(R? - R2) 



TO2e2(f-A) 

R?Ri 

^2g2($-A) 

R?Ri 



(g2-l)(g2^26)2 

(Ri+R2)2 

(l^g^)(z' + 26)^ 
(Ri+R2)2 



(22c) 
(22d) 
(22e) 



where k = ak, Ri = aRi and R2 = aR2. 
From H22a|) and the condition 

^ ^ _fc(R[+R|^ ^ ^ 
(Ri - R2)R?Ri " ' 

we have p > 0, i.e., the strong energy condition is satis- 
fied. Also from R^ - R| = Ak{h + &), with /i(i) = ^2/2, 
we can show that 



A = 



< 



(Ri + R2)(Rf + Ri) 
4(/i + 6)R2r2 

R? Ri + 4fc(/i + 6) 



{h + 6)R2 (/i + 6)r2 
Also h + b> b and R2 > ^ — give us 
, (b-k)^+4kb 

&(&- fc)2 

Then the condition A < 1 can be cast as 

{b + k f <b{b-k)\ (23) 

that leads also to b ^ k. This last condition assure the 
nonsingular behavior of the energy density e and the az- 
imuthal stress p^. 

From Ij22b|) and ?- = we have that e > 0. When z ~ 1 
the condition 

^_ 2mf2(R^+R^) ^ ^ 

R2R2[(Ri+R2)2_4fc2] - ' 

gives us e . Since Ri > R2, Ri + R2 > 25, R2 > r 
and Ri > 6 + A: we have, 

m 

B 



Then the condition B <\ leads to 

0<m< (62-P)(6 + fc). 



(24) 



This relation yields also 6 > A: as a condition to have 
TO > 0. 

For any other value of r and z is not easy to obtain 
constraints over the parameters fc, b and to in order to 
have e > 0. The analysis is better done graphically. By 
considering different values of k and b that fulfill the con- 
dition l|23|) we find that e > everywhere in the disks only 
if we take for fh a value less than a tenth of the upper 
limit provided by the condition (|24|l . As an example, in 
Fig. Elwe plot the dimensionless densities p and e for a 
disk with to = 3, 5 = 3.5 and k = 1. We can see that, as 
in the Chazy-Curzon disk, the density have a maximum 
at the center of the z = plane and then it decreases 
monotonously as r increase and also for z ±a. Also p 
and e have similar magnitudes. 

The behavior of the stresses is also similar to the pre- 
sented in the Chazy-Curzon disk. Again, is better to do 
a graphical presentation. In Fig. \7\ we plot the dimen- 
sionless azimuthal stress and the vertical pressure Pz 
for the disk with m = 3, 6 = 3.5 and k = 1. Again we 
have that pr = —pz- As with the Chazy-Curzon disk, the 
azimuthal stress is negative at the central region of the 
disk, then increase to have a positive maximum at the 
z = plane, for a value of r « 2.5a, and finally decrease 
monotonously. The behavior of the vertical pressure is 
like of the densities, with a maximum at the center of 
the z — plane and then monotonously decreasing for 
increasing r. Also, p^ = for z — ±a. From Fig. Eland 
[3 we have 



Pv 
e 



{b^-k^){b + k) 





Pr 




Pz_ 




e 




e 


and so the disks are 


also 


in 



< 0.2, 



< 0.1, 
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FIG. 6: For a thick disk obtained from the Zipoy-Voorhees solution with m = 3, 6 = 3.5 and fc = 1 we plot (a) the effective 
Newtonian density p and (b) the energy density I. 





FIG. 7: For a thick disk obtained from the Zipoy-Voorhees solution with in = ?>, b = 3.5 and fe = 1 we plot (a) the azimuthal 
stress ptp and (b) the vertical pressure p^. 



nant energy condition. Thin disks based in the Zipoy- 
Voorhees metric were considered in Ref. H. 



IV. THICK DISKS FROM DE 
SCHWARZSCHILD METRIC IN ISOTROPIC 
COORDINATES 

For a static spherically symmetric spacetime the metric 
in isotropic spherical coordinates (i, i?, 9, ip) can be cast 
as 

ds^ = -e^'^dt^ + e^^[dR^ + R^{d9^ + sin^ Odif^)], (25) 

where $ and A are functions of R only. In isotropic 
cylindrical coordinates (t, ip, r, z) the metric (I25|) takes 
the form 



dinates, 



ds^ = 



^^•^dt^ + e^^ir^dp^ + dr' + dz% 



(26) 



where now $ and A depends on r and z. 

We will now apply the "displace, cut, fill and reflect" 
method to the Schwarzschild solution in isotropic coor- 



In 



2R- m 



A = In 



2R + m 



1+2R 



(27a) 



(27b) 



with TO a positive constant and = -I- z^. Now we 



put R^ 



{h + 5)^ where 6 is a positive constant and 



h{z) is given by Q). 

From the Einstein equations and the orthonormal 
tetrad {\/", X", F", Z''}, where 



V 
X" 



(1,0,0,0) 



— (0,1,0,0) 
r 



'(0,0,1,0) 



no, 0,0,1) 



(28a) 

(28b) 
(28c) 
(28d) 
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we find that the energy-momentum tensor of the disk can 
be written as 

Tab - eVaVb + p^XaXb + PrYaYb + p.ZaZb . (29) 

Here e — — is the energy density, p^p = is the 
azimuthal stress, that is equal to the radial stress pr = 

I 



T^, and Pz = is the vertical stress. The effective 
Newtonian density is given by p = e + 2pj. + = e + 
2p<^ +Pz- 

From (|27|l we obtain, using the dimensionless variables 
previously defined. 



Pv 



Pr 



3mR 



2(3F + 2{h - 1))R2 + i{P + 26)2(1 - z^) 



3m 



(2R-m)(2R + m)5 
2(3z2 + 2(6 - 1))R2 + 3(^2 + 26)2(1 - p) 



16m2 



2(2R + m)5 
(3^2 + 2(6 - 1))R2 + (i2 + 26)2(1 - 



16m 



(2R- m)(2R + 7fi)5 
(3i2 + 2(6 - 1))R2 + (52 + 26)2(1 - 



(2R- m)(2R + m) 
16m2(f2 + 26)2(1 -52) 
(2R-?fi)(2R + m)S ' 



(30a) 

(30b) 

(30c) 

(30d) 
(30e) 



From (|30b|l and 6 > 1 we have e > 0. On the other 
hand, from H30a|) and 2R > m we have p > 0. Since 
R > 6, this last condition is equivalent to 26 > rh. There- 
fore, when 6 > 1 and < to < 26 we will have disks in 
agreement with the weak and strong energy conditions. 
Also this values of to assure the nonsingular behavior of 
p, Pip, Pr and Pz- We also have p^ — pr > and Pz > 0. 
The vertical and horizontal stress are then pressures. 

As in the previous section, we perform a graphical anal- 
ysis of the solution. In Fig. |S|we plot p and e for a thick 
disk obtained from the Schwarzchild isotropic solution 
with fh — h = 1. The horizontal and vertical pressures 
are plotted in Fig|51 All the four quantities have a simi- 
lar behavior, with a maximum at the center of the z = 
plane and then monotonously decreasing with increasing 
r and z. The relative magnitudes of the densities and 
pressures are such that p > e > Pp = Pr ~ Pz- We have 







Pr 




El 


e 




e 




e 



and so the disks are in agreement with all the energy 
conditions. Thin disks based on the Schwarzschild solu- 
tion in isotropic coordinates were studied in 29 ] , whereas 
thin disks based in the Schwarzschild metric in Weyl co- 
ordinates were studied in For the disks in isotropic 
coordinates we have matter with radial pressure equal to 
the azimuthal pressure (isotropic matter) an for the disks 
in Weyl coordinates we have zero radial pressure. 



V. DISCUSSION 

We presented a method to obtain exact general rela- 
tivistic thick disks as a generalization of the "displace, 
cut and reflect" method commonly used to obtain New- 
tonian and relativistic thin disks. The generalization was 
be done by means of the transformation z — > h(z) + 6, 
where h(z) is an even function of z and 6 is a positive 
constant. The function h{z) must be selected in such a 
way that the metric tensor and and its first derivatives 
will be continuous across the plane z = 0. 

All the cases considered leads to thick disks with simi- 
lar behavior of the energy and Newtonian effective densi- 
ties: a maximum at the center of the central plane of the 
disks, the z — plane, and then monotonously decreasing 
for increasing r and z. 

We found that, when the method is applied to vacuum 
Weyl spacetimes, the thick disks presents radial tension 
and vertical pressure. The azimuthal stress is negative 
at the central region of the disks, then have a positive 
maximum and finally decrease monotonously for large 
values of r and z — ±a, where 2a is the thickness of the 
disk. The disks obtained are in full agreement with all 
the energy conditions. 

On the other hand, when the method is applied to the 
Schwarzschild isotropic metric all the stresses are pres- 
sures and have a behavior like the densities. The disks 
obtained are also in full agreement with all the energy 
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FIG. 8: For a thick disk obtained from the Schwarzschild isotropic solution with m = 6 = 1, we plot (a) the effective Newtonian 
density p and (b) the energy density e. 




FIG. 9: For a thick disk obtained from the Schwarzschild isotropic solution with m = 6 = 1, we plot (a) the horizontal pressure 
Pip = pr and (b) the vertical pressure pz- 



conditions. 

We plan to extend the models of thick disk presented 
along this lines by considering more elaborate functions 
h{z) and by the incorporation of new properties like ro- 
tation, either electric or magnetic fields or both. Also 
we believe that the study of stability in this disks can 
produce some non trivial results. 
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